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ABSTRACT 

We have detected narrow Hi 21cm and Cl absorption at z L4- 1.6 towards Q0458-020 and Q2337-011, 
and use these lines to test for possible changes in the fine structure constant a, the proton-electron mass ratio 
^, and the proton gyromagnetic ratio gp. A comparison between the Hi 21cm and Cl line redshifts yields 
AX /X = [+6.8 ± 1 .0] X 10"^ over < (z) < 1 .46, where X = gpa^/^j,, and the errors are purely statistical, from 
the gaussian fits. The simple line profiles and the high sensitivity of the spectra imply that statistical errors 
in this comparison are an order of magnitude lower than in previous studies. Further, the Cl lines arise in 
cold neutral gas that also gives rise to Hi 21cm absorption, and both background quasars are core-dominated, 
reducing the likelihood of systematic errors due to local velocity offsets between the hyperfine and resonance 
lines. The dominant source of systematic error lies in the absolute wavelength calibration of the optical spectra, 
which appears uncertain to ~ 2 km/s, yielding a maximum error in AX/X of ^ 6.7 x 10"^. Including this, we 
obtain AX/X = [+6.8 ± I. Q(statistical)±6.1{max. systematic)] x 10"^ over < (z) < 1.46. Using literature 
constraints on A/i/ /i, this is inconsistent with claims of a smaller value of a from the many-multiplet method, 
unless fractional changes in gp are larger than those in a and /i. 

Subject headings: atomic processes — galaxies: high-redshift — quasars: absorption lines 



1. INTRODUCTION 

A critical assumption in the standard model of particle 
physics is that low-energy coupling constants and particle 
masses do not vary with space or time. This assumption 
breaks down in most theories that atte mpt to uni fy the stan- 
dard model and general relativity (e.g. Marcianol[l984) . The 
detection of such spatio-temporal variation in coupling con- 
stants like the fine structure constant a, or the ratios of parti- 
cle masses (e.g. the proton-electron mass ratio /i = nip /me), 
would imply new physics beyond the s tandard model and is 
hence of great interest (e.g. IUzanll2003h . 

Comparisons between the redshifts of spectral lines de- 
tected in distant galaxies provide an important tool to probe 
changes in a, /i and the proton gyromagnetic ratio g,, over 
cosmological ti mes (e. g. [Wolfe et al. 1976; Dzuba et al. 1999t 
IChengalur & Kanekad 12003; Flambaum & Ko/lov 200J. 
Most such studies have yielded constraints o n changes i n a, ij , 
and gp, with different systematic effects [see iKanekaii (l2008h 
for a recent review]. At present, the only technique that has 
found statistically-significant evidence for a change in one 
of the fundamental constants is the "many-multiplet" method 
dPzuba et al.ll999h . iMurphv et al.1 ( 120041) obtained [Aa/a] = 
(-5.7 ± 1.1) X 10"'' from Keck High Resolution Ultraviolet 
Echelle Spectrometer (HIRES) optical spectra of 143 ab- 
sorbers with a mean redshift (z) = 1 .75, suggesting that a was 
smaller at earlier times. Other studies, applying a similar tech- 
nique to Very Large Telescope (VLT) data on smaller samples , 
have not confirmed this result (e.g. 'Levshakov et al.l 120061: 
ISria nand et al. 2007). However, Murphy et al. (2008b) argue 
that the errors in these studies have been under-estimated, and 
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that results from small samples are more prone to systematic 
effects related to the fitting of spectral co mponents to com- 
plex absorption profiles. In the case of n, iKing et al.1 (l2008l) 
used VLT Ultraviolet Echelle Spectrograph (UVES) spectra 
of ro-vibrational H2 lines in three damped Lya systems at 
z ~ 2.6-3 to find [A^/m] = (-2.6 ±3.0) x 10-^ Note that 
none of the above error estimates include recently-detected 
systematic effects due to distortions in the wavel ength scales 
of the HIRES and UVES spectrographs (e.g. iGriest etaP 
l201(]h . Finally, constraints on changes in a, fi and gp have 



also been o btained f rom radio techniques (e.g. iCarilli et al 



'2000; KanekareLalJL2004; Kan ekar et al.ll2005l;lMurphv et al 
l2008a), although at lower redshifts. 

Comparisons between the redshifts of Hi 21cm (hyperfine) 
and ultraviolet resonance dipole transitions are sensitive to 
changes in X = gpO^/fi (Wolfe et al. 1976). The best res- 
onance lines for this method are those arising from neutral 
atomic species (e.g. Cl, Fel, etc), as these species are most 
likely to be physically associated with the Hi. The Cl mul- 
tiplets are likely to be the best among the neutral resonance 
transitions as they typically arise in cold gas which also gives 
rise to the H i 21cm absorption (e.g. Jenkins & Tripp 200ft 
ISrianand"eral. 2005) . The ionization potentials of Cl and Hi 
are also similar, 1 1.3 eV for Cl and 13.6 eV for Hi (cf. Mgl, 
which is more easily detectable than Cl in high-z absorbers, 
has an ionization potential of 7.6 eV, as well as a high dielec- 
tronic recombination rate that can give significant Mgl ab- 
sorption in warm, ionized gas; Pettini et al. 1977). Finally, 
absorbers with a single (or dominant) spectral component in 
both Hi 21cm and Cl transitions are best-suited for such stud- 
ies. 

Hi 21cm and Cl absorption have hitherto both been de- 
tected in only one high-z absorber, the z ^ 1.776 system 
towards 1331+17 0, yielding a weak cons traint on changes 
in X = gpa^/^JL (ICowie & Songailal [l995h . However, the 
Cl line in this absorber ha^s two clear spectral components 
(iDessauges-Z avadsk y et al.ll200 4'). implying ambiguities (and 
thus, large systematic errors) in the comparison with the 
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Hi 21cm line. In this Letter, we report the detection of nar- 
row, single-component Hi 21cm and Cl absorption in two ab- 
sorbers at z 1.4-1.6, that allow a high-sensitivity study of 
changes in the fundamental constants. 

2. SPECTRA AND RESULTS 

The Giant Metrewave Radio Telescope (GMRT) and the 
Green Bank Telescope (GBT) were used to detect Hi 21cm 
absorption at z ~ 1 .3609 towards Q2337-01 1 and z 1 .5605 
towards Q0458-020, respectively, in a survey for Hi 21cm ab- 
sorption in strong Mgll absorbers [GMRT proposals 7NKa02, 
10NKa02, and GBT proposal 6A-026; see Kanekar etal. 
(|2009) for details]. The Hi 21cm spectra are shown in the 
lower panels of Figure[T] The root-mean-square (RMS) noise 
values, measured from off-line regions in the spectra, are 
- 0.03 per 1.9 km/s channel (Q2337-01 1) and - 0.0037 per 
1.3 km/s channel (Q0458-020), in optical depth units. 

We then carried out a search for redshifted Cl absorp- 
tion from the two absorbers, using the high-sensitivity 
Keck-HIRES spectrum of Prochaska & Wolfe ( 1997) towards 
Q0458-020 (observed on 31 October and 1 November 1995), 
and a new HIRES spectrum of Q2337-011 (observed on 18 
and 19 September 2006). This resulted in the detection of 
the ClA1560 and ClA1657 ti-ansitions at z 1.3609 towai'ds 
Q2337-011, and the ClA1560 and Cl*A1561 transitions at 
z - 1.5605 towards Q0458-020. Of these, the ClA1560 Hne 
towards Q2337-011 is blended with another line, while the 
Cl*A1561 hne towards Q0458-020 is detected at low sig- 
nificance; these will hence not be used in the later analysis. 
The CIA1657 Une towards Q2337-011 and the ClA1560 Une 
towards Q0458-020 are shown in the upper panels of Fig- 
ure [T] The signal-to-noise ratios per pixel are 25 at a res- 
olution of - 50000 (Q2337-01 1), and - 15 at - 37000 
(Q0458-020), in the vicinity of the above Cl transitions. 

For all spectra, the RMS noise was measured from 
absorption-free spectral regions around the line in question. A 
single-gaussian model was then used to independently fit each 
of the 1-D spectra in the Hi 21cm and Cl lines. This model 
yielded an excellent fit to all spectra, with reduced chi-square 
values < 1 .06 in all cases, and no evidence for statistically- 
significant features in the residual spectra. A Kolmogorov- 
Smirnov rank-1 test found all residual 1-D spectra to be con- 
sistent (within 1.1 fT significance) with being drawn from a 
normal distribution. The RMS noise values used for the spec- 
tral fits were scaled (marginally) to obtain = 1 in all cases; 
the resulting fits were then used to measure the peak absorp- 
tion redshift for each transition. These redshifts are listed 
in columns (2) and (3) of Table [U note that the ClA1560 
and ClA1657 lines have laboratory vacuum wavelengths of 
1560.3092 A and 1657.9283 A, respectively (lMortonll2003"h . 
while t he Hi 21cm Un e frequency is 1420.405751766 (1) 
MHz (Essen et al. 1971). 

A useful test of the limiting accuracy in such fits can be ob- 
tained by combining the minimum uncertainties contributed 
by all line pixels to determine the best velocity accur acy that 
might be obtained with a given spectral fit (Murp hy et al.l 
I2008bi) . In all cases, the redshift errors in Table [1] are larger 
than this limiting accuracy (by factors < 2), as expected for 
real spectra. We also tested that adding additional compo- 
nents to the fits does not significantly alter the results. 

Assuming that the Hi 21cm and Cl lines arise in the 
same gas, the fractional change in A" = gpO^ / ji is related to 
the observed Hi 21cm and Ci redshifts by AX/X = [zqi- 



TABLE 1 

Parameters of the single-gaussian fits to the Hi 21cm and Ci 
spectra towards q2337-()l 1 and q0458-020. 





Z21cm 


^Cl 




AX/X 






km/s 




Q2337-011 


1.3608644 (13) 


1.3608801 (15) 


1.99 ±0.25 


+6.64 ±0.84 


Q045 8-020 


1.5605300 (25) 


1.5605480 (38) 


2.11 ±0.53 


+7.0 ± 1.8 



t The velocity off.set of the Cl line from the Hi 21cm line. 

Z2icmV[l +3, where z= (zQ-Z2icm)/2 (I^EiEillllTll). 

This yields AX/X = (+6.64 ± 0.84) x 10"^ for the z - 1 .3609 
absorber towards Q2337-011 and AX /X = (+7.0 ± 1.8) x 
10"'' for the z 1.5605 absorber towards Q0458-020. Av- 
eraging these values (with equal weights) gives the result 
AX/X = (+6.8 ± 1.0) X lO"*", over < z < 1.46. We em- 
phasize that the errors quoted here are purely statistical ones, 
from the fits; systematic effects are discussed in the next sec- 
tion. 

3. SYSTEMATIC EFFECTS 

At the outset, it should be emphasized that, unlike the 
many-multiplet method which only requires accurate relative 
wavelength calibration between different lines in the same op- 
tical spectrum, the comparison between Hi and Cl lines re- 
quires accurate absolute wavelength calibration of the opti- 
cal spectra. On the other hand, the many-multiplet method 
is based on a first-order effect (the wavelengths of transitions 
used in this analysis have the same zeroth-order dependence 
on a), while the present method uses the zeroth-order depen- 
dences of the hyperfine and resonance line frequencies on a, 
1^1 and gp. This implies that systematic effects (e.g. due to 
wavelength mis-calibration) are less important by an order of 
magnitude in the hyperfine/resonance comparison than in the 
many-multiplet analysis. For example, a velocity uncertainty 
of 0.3 km/s implies an error of [Aa/a] = 10"^ in the many- 
multiplet method dMurphv et alj|2001h . but of AX/X = IQ-^ 
in the hyperfine/resonance comparison. 

The Cl lines in both absorbers are at higher redshifts than 
the Hi 21cm lines, by ^ 2 km/s. Systematic effects that might 
cause such a velocity offset include (1) different relative iso- 
topic abundances at z ~ 1.5 from the local Universe, that 
might cause the rest Cl wavelengths in the two absorbers to be 
higher than their laboratory values, (2) errors in the absolute 
wavelength calibration of the optical spectra, and (3) "local" 
velocity offsets between the Cl and Hi 21cm lines in the two 
absorbers. Note that frequency mis-calibration of the Hi 21cm 
spectra is unlikely to be the source of such errors, as the abso- 
lute frequency scale here is set by the accuracy of masers and 
local oscillators (typically < 10 Hz, two orders of magnitude 
lower than the observed velocity offset between the Hi 21cm 
and Cl lines). 

The velocity shifts of the '"^C and '^'C isotopic transi- 
tions relative to the main '^C isotopic lines are, respec- 
tively, +0.46 km/s and +0.84 km/s for the ClA1657 line, 
and -3.29 kni/s and -6.10 km/s for the ClA1560 line 
dBerengut et al.1 12006'). The isotopic velocity shift is clearly 
too small to account for the observed velocity offset in the 
case of the ClA1657 line. Conversely, in the case of the 
ClA1560 line, the isotopic shift yields the opposite sign from 
the observed velocity offset between Cl and Hi lines. We 
can thus rule out the hypothesis that differing carbon isotopic 
abundances (compared to Galactic values) might account for 
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Fig. 1. — Hi 21cm and Cl spectra towards [A] Q2337— Oil and [B] Q0458-020. In each panel, the solid line shows the 1-gaussian fit to the spectrum. The 
small-dashed and large-dashed vertical fines indicate the Cl and Hi 21cm redshifts, respectively; the Cl redshift is seen to be higher than the Hi 21cm redshift in 
both panels. 

duce errors in the wavelength scale (Murphy et al. 2001). Fur- 
ther, while the data towards 02337-01 1 were obta i ned in 
2007, with the use of the image rotator, iGriest etaTI (120101) 
have found evidence for drifts in the Keck-HIRES wavelength 
scale with time, with amplitudes of ^ 2 km/s over multiple 
observing epochs. The source of these velocity drifts is still 
unclear, although temperature changes, changes in the posi- 
tion of the quasar on the slit, as well as physical shifts in the 
echelle grating or cross-disperser are all possible causes. 

To test for systematic effects in the absolute wavelength 
calibration of the optical spectra, we analysed archival Keck- 
HIRES data towards 0458-020 from 5 and 6 October 2004 
(obtained using the image rotator). While this spectrum 
has lower S/N than the 1995 spectrum (as the optically- 
variable quasar was in a fainter state in 2004), the unsaturated 
NillA1317 ti-ansition from a higher-redshift (z = 2.03945) 
damped Lyman-a system is clearly detected in both spec- 
tra; this Une lies at 4003 A, within 10 A of the ClA1560 
line shown in Figure[TtB]. Cross-correlating the Nill lines de- 
tected in the two HIRES spectra yielded a cross-correlation 
peak at a velocity offset of +1.60 ±0.84 km/s; the error 
was determined by cross-correlating 10000 pairs of simulated 
spectra with the shape and noise properties of the actual spec- 
tra. If this velocity offset between the Nin lines arises due 
to atmospheric dispersion across the slit during the 1995 ob- 
servations, it would imply that the Cl line redshift has been 
under-estimated in the 1995 spectrum. Correcting for this ef- 
fect would increase the redshift offset between the Hi 21cm 
and Cl lines. 

We also separately analysed the Keck-HIRES data to- 
wards 02337-011 from 18 and 19 September 2006, to test 
whether the same line redshift was obtained from the two 
runs. Figure |2] shows the ClA1657 line profiles obtained 
on 18 and 19 September; these are clearly offset from each 
other. A single-component gaussian fit was used to measure 
the ClA1657 redshift from each spectrum; the offset between 
the ClA1657 redshifts is 2.15 ±0.44 km/s. Similar velocity 
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Fig. 2. — A comparison between the CIA1657 line profiles obtained on 18 
and 19 September, plotted as squares and triangles, respectively. The dashed 
vertical lines indicate the Cl redshifts obtained from 1-gaussian fits to each 
spectrum; these are offset from each other by ~ 2 km/s. The solid circle 
shows the Hi 21cm redshift; see main text for discussion. 

the observed velocity offset between the Hi 2 1 cm and Cl tran- 
sitions. 

Mis-calibration of the absolute wavelength scale of the 
optical spectra coul d arise from a number of causes (e.g. 
iMurphv et al.l 120011) . In particular, the spectrum towards 
00458-020 was obtained in 1995, before Keck-HIRES was 
fitted with an image rotator. It was thus not possible to hold 
the slit perpendicular to the horizon during these observations, 
implying that atmospheric dispersion across the slit could pro- 
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offsets (^2 km/s) were seen between other lines in the two 
spectra. Note that a single ThAr lamp exposure was used to 
calibrate the two science exposures on Q2337-011 on each 
day; the ThAr exposure was taken immediately before/after 
the science exposures. Interestingly, the ClA1657 redshift of 
18 September [zqi = 1.3608694(26)] is in reasonable agree- 
ment with the Hi 21cm redshift [z2Lcm = 1-3608644(13)] , 
indicated by the solid circle in Figure[2l the redshift offset be- 
tween the Hi 21cm and Cl lines is dominated by the Cl data 
from 19 September The offsets of ^ 2 km/s between our two 
observing epochs are consistent wit h the velo city drifts in the 
HIRES wavelength scale found by iGriest et ai.l (2010) (see 
their Figure 5). We hence conclude that the absolute wave- 
length scale of HIRES could be in error by ~ 2 km/s. 

Finally, "local" velocity offsets between the Hi 21cm and 
Cl lines within the absorbing galaxies might also contribute 
to differences in the measured line redshifts. Both Cl and 
Hi 21cm absorption are expected to arise in cold gas, and, 
for both background quasars, a significant fraction of the ra- 
dio flux densi ty at the Hi 2 1 cm lin e frequency arises from a 
compact core (iKanekar et al.l2069h . In fact, Q2337-01 1 has a 
highly inverted spectrum, indicating that most of the flux den- 
sity arises from a self-absorbed core (Kanekar et al. 2009). 
The cold cloud producing Cl absorption against the optical 
quasar is thus also likely to give rise to Hi 21cm absorption 
against the radio core. Assuming that gravity is importan t 
for cloud confinement on small scales. Mur phy et al.l (l2003h 
obtain a velocity dispersion of 0.1 km/s between different 
species in an individual cloud, which would yield a systematic 
error of ~ 3.3 x 10"^ in AX/X, significantly lower than our 
statistical errors. We note, however, that larger velocity off- 
sets due to small-scale structure in the absorbing gas cannot 
formally be ruled out. 

4. DISCUSSION 

Prior to this work, the most sensitive result constraining 
changes in X = ^p O^/iJ, using the hyperfine/resonance com- 
parison was that of iTzanavaris et al. (2007). These authors 
compared redshifts of the deepest absorption in Hi 21cm and 
low-ionization metal lines to obtain AX/X = (6.3 ±9.9) x 
10"^ from a sample of nine absorbers at 0.23 < z < 2.35. 
Note that these error estimate do not include systematic ef- 
fects. Considering only statistical errors, our result, AX/X = 
[+6.8 ±1.0(statistical)] x 10"^ over < (z) < 1 .46, is an or- 
der of magnitude more sensitive than that of ITzanavaris et all 
(120071). Furthe r , the low-ionization metal lines used by 
ITzanavaris et al.l (120071) could also arise in warm Hi or ion- 
ized gas, and are not necessarily associated with the cold Hi 
that gives rise t o the Hi 21cm ab sorption. Most of the ab- 
sorbers used by ITzanavaris et al.l ([2007) also have complex 
Hi 21cm and metal line profiles, and it is not necessary that 
the deepest absorption in the two types of transitions arises 
in the same spectral component (Kanekar et al. 2006). This 
could give systematic err ors of > 10 km/s, far larger than the 
statistical errors of Tzana varis et al.l (120071) . or the systematic 
errors in the present result. 

The comparison between hyperfine and resonance transi- 
tions directly probes changes in X = gpa^/ii, and one can- 
not obtain independent constraints on the individual con- 
stants without additional assumptions. Our result gives 2 x 
[Aa/a] + [Agp/gp] - [A/i//i] = [+6.8 ± l.Q{statistical) ± 
6.1 {systematic)] x 10"^. Note that 6.7 x 10"^ is the maxi- 
mum estimated error due to systematics in the absolute wave- 



length calibration of the optical spectra, (and not a la es- 
timate, as in the case of the statistical error). This implies 
that 2 x [Aa/a]+[Agp/g p]- [A/z//z] > [+0.1 ±1.0] x 10"^ 
The Keck/HIRES result of lMurphy et aPdlOOl is [Aa/a] = 
(-5.7 ± 1.1) X 10"^, with (z) = 1.75. Consistency between 
the two results would require either (1) an additional wave- 
length calibration error of ^ 1.7 km/s in the Cl spectra, which 
appears unlikely, (2) a "local" velocity offset of ^ 1 .7 km/s 
(and with the appropriate sign) between the Hi and Cl lines 
in both absorbers, (3) under-estimated errors in the many- 
multiplet result [e.g. the distortions in the wavelength scale 
found by IGriest etaH (120101) 1. or (4) [Agp/gp] - [A^i/ ^j] > 
(+1.14 ±0.24) X 10"^ at (z) = 1.46. In other words, assum- 
ing that the errors in one (or both) of the results have not been 
under-estimated, consistency between the results requires that 
fractional changes in and/or gp are comparable to those in 
a. Strong constraints are available on fractional changes in 
/i at both higher and lower redshifts, [A/i/ /i] < 1 .6 x 10"^ at 
z - 0.685 (Murphy et al."2008a') and [A^/fj] < 6.0 x 10"'' at 
z 2.8 (King et al. 2008), making it unlikely that [A^/^] 
10"^ at z ~ 1.46. We thus conclude that the present result ap- 
pears inconsistent with the smaller value of a at (z) ~ 1.75 
found by Murphv et al. (2004), unless fractional changes in 
the proton gyromagnetic ratio gp are larger than those in a 
and /i. 

In summary, we have detected narrow Hi 21cm and Cl ab- 
sorption in two absorbers at z 1 .4- 1 .6 towards Q0458-020 
and Q2337-011, using the Keck telescope, the GMRT and 
the GBT. The Cl and Hi 21cm line frequencies have differ- 
ent dependences on the fundamental constants a, fi = nip/m^ 
and gp, allowing us to use the measured Cl and Hi 21cm 
line redshifts to test for putative changes in these constants. 
Comparing the Hi 21cm and Cl redshifts in the two ab- 
sorbers yields the result AX/X = [+6.8 ± l.Q(statistical)± 
6. 7 (systematic)] x 10~*'overO < (z)1.46, where X =gpQ^ / fi. 
This is inconsistent with evidence for a smaller value of a 
at similar redshifts from the many-multiplet method, unless 
fractional changes in gp are larger than those in a and fi. Sys- 
tematic errors in the hyperfine/resonance comparison are cur- 
rently dominated by errors 2 km/s) in the absolute wave- 
length calibration of the optical spectra. However, the com- 
parison between Hi 21cm and Cl lines has a high sensitivity, 
and it should be possible to significantly reduce systematic ef- 
fects by the us e of new cal ibration techniques (e.g. laser fre- 
quency combs: ISteinmetz e t al. 2008). Increasing the number 
of detections of redshifted, narrow Hi 21cm and Cl absorption 
is thus of much importance. 
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